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Variation of the 77 opacity by the He 11 Lyman continuum constrains 
the location of the 7-ray emission region in the blazar 3C 454.3 
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ABSTRACT 

We study spectral properties of the brightest 7-ray blazar 3C 454.3 using 138 weeks of ob- 
servations by the Fermi Gamma-ray Space Telescope (Fermi). We probe the behaviour of the 
source as a function of time at different brightness levels. The Fermi spectra in the GeV range 
can be well described by a wide underlying lognormal distribution with the photon-photon ab- 
sorption breaks produced by the He II and H I Lyman recombination continua (LyC). We find 
a power-law dependence of the peak energy on flux and discover anti-correlation between the 
column density of the He II LyC and flux. This implies that the 7-ray emission zone lies close 
to the boundary of the high-ionization part of the broad-line region and moves away from the 
black hole when the flux increases. Identification of the 7-ray production with the relativistic 
jet, implies that the jet is already accelerated at sub-parsec distances from the central black 
hole, which favours the Blandford-Znajek process as the jet launching mechanism. 

Key words: black hole physics - galaxies: active - galaxies: jets - gamma rays: observations 
- quasars: emission lines - quasars: individual (3C 454.3) 



1 INTRODUCTION 

The flat spectrum radio quasar (FSRQ) 3C 454.3 at redshift z — 
0.859 is by far the brightest blazar in the 7-ray range during the life- 
time of the Fermi Gamma-ray Space Telescope (Fermi). Since the 
start of operation in August 2008 the Fermi Large Area Telescope 
(LAT) has detected more than hundred thousand photons from this 
object. Its brightness has been changing by at least two orders of 
magnitude (see Fig.fJJ: after a bright state in 2008, the blazar turned 
to a very low state in the spring of 2009, then the flux has been grad- 
ually increasing. The bright flares then occurred in November 2009 
and April 2010 and the brightest 7-ray flare ever observed from a 
blazar happened in November 2010. 

The early data from the first 160 days of the Fermi observa- 
tions have demons trated a clear spectral break around 2.5-3 GeV 
dAbdoetal.ll2009h . Similar breaks in the GeV ranee have been 
observed in several other FS RQs and low-sync hrotron peak BL 
Lac o bjects 1 Abdo et all 2010li . Later works by lAckermann et al.l 
d2010h and I Abdo et alj J2OI II) dedicated to 3C 454.3 found the 
spectral hardness-flux correlation at a nearly constant break energy. 
First inteipretations of the breaks invoking the external Compton 
mechanism with a trun cated electron spectrum (Abd o et al1l2009l : 
iFinke & Dermeij|2010h do not reproduce as sharp a break as ob- 
served. The stability of the break energy favours instead the absorp- 
tion of the GeV photons by the photon-photon (77) pair production 



E-mail: stern.boris@gmail.com, juri.poutanen@oulu.fi 



on the He II and H I Lyman recombinatio n continua (LyC) from the 
broad-line region (BLR), as proposed bv lPoutanen & Sternl ( l2010l 
hereafter PS 10). A relatively high opacity in the He II LyC observed 
in a number of bright blazars implies the location of the 7-ray emit- 
ting region within the highly ionized inner part of the BLR. 

Motivated by the extremely rich photon statistics from 
3C 454.3, in the present work we particularly try to answer the 
following questions: 

(i) Is the spectral break consistent with the 77-absorption on the 
He II LyC at various source fluxes? We do not have a freedom of 
choosing the position of the break, with the only free parameter 
being the LyC column density. A good fit at a rich statistics would 
be a strong argument in favour of the absorption hypothesis. 

(ii) Is there any correlation between the break sharpness (i.e. ab- 
sorption opacity in our interpretation) and the flux? If yes, then this 
would imply a varying location of the 7-ray emitting region. 

(iii) What is a proper characterization of the spectrum below 
the break? In PS 10 the hypothesis for the underlying (unabsorbed) 
spectrum was a power-law. The underlying spectrum can actually 
be slightly curved with respect to the power-law. 



2 DATA AND THEIR ANALYSIS 

The light curve of 3C 454.3 from the start of the Fermi data records 
to 2011 April 1 is shown in Fig.fT] We split the data into ten time 
intervals covering periods of more or less constant flux levels and 
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Figure 1. The light curve of the 7-ray blazar 3C 454.3 from the beginning of the Fermi data records to 2011 April 1. Horizontal bars mark the time intervals 
used in the spectral analysis (see Table[T}. 



analyze the spectrum in each interval. We use our own software 
for the data analysis because the standard software distributed by 
the Fermi LAT team does not contain the required spectral models. 
The energy-dependent exposure function was calculated using the 
spacecraft pointing history: 



Exposure(.E, AT) 



S(E,0(t)) dt, 



(1) 



where AT is the time interval of interest, 6 is the angle between the 
detector axis and the direction to the object and S(E, 8) is the de- 
tector effective area at energy E. The spectral models folded with 
the exposure function are fitted to the observed counts in 13 loga- 
rithmically spaced energy bins in the energy range 0.15-60 GeV. 

We included the events of classes 3 and 4 and imposed the 
cuts at zenith angle <105°. We used P6_V3_DIFFUSE ver- 
sion of the response function. The background was measured in 
a circle with the radius 6°, separated from the object by ~15° 
and avoiding local sources. We accumulate counts in the circle 
centered at the source location with the energy-dependent radius 
r = max{10°(T/100 MeV) -0 ' 7 , 0?6} (which corresponds to 
95% containment). Statistical errors were treated as Gaussian, ex- 
cept in a few bins at higher energies, where the number of photons 
is low, we use Poisson likelihood adding —2 In P(n, /x) to y 2 (here 
n is the number of counts in the bin and fi is the prediction of the 
model). The number of such bins is small and the meaning of \ 2 is 
not significantly affected. For the minimization we use the standard 
code MINUIT from the CERN library. We checked that our software 
reproduces within statistical errors the results of the gtlike task for 
standard spectral models such as power -law and broken power-law 
(compare table 2 in PS 10 with table 1 in lAbdo et al.l2010l) . Fig.[Ha) 
presents the spectra of 3C 454.3 accumulated in each time interval. 



3 SPECTRAL FITS 

First, we fit the data with a simple power-law model. This is our 
first null-hypothesis against which we check the significance of the 



break. The resulting values of \ 2 are given in column 3 of TableQ] 
Then we add the 77-absorption by photons in two "lines" at 54.4 
and 13.6 eV (double absorber model, DA), which correspond to 
the energies of the He II and H LyC, which are the strongest fea- 
tures in the BLR spectra at relevant energies. These continua are 
basically line-like because the photoionized region temperature is 
typically much below the corresponding ionization potentials of H 
and He. Within this model, the spectral breaks at expected at the 
energies 261/(Tii nc (eV)[l + z]) GeV corresponding to the thresh- 
old of the 77 absorption on photons of energy £qi nc and redshifted. 
In 3C 454.3, the breaks should a ppear at 2.6 and 10.3 GeV, which 
is con sistent with observations dAbdo et al.l2009l ; lAckermann et al.l 
l20ld :PS10; seeFig.H. 

The absorption strength by the isotropic line photons is pa- 
rameterized by the Thomson optical depth (see PS 10): 

Lllne,45 10 eV 



TT 



N ph a T « 110 



(2) 



Rl8 £"line 

where N p ^ is the photon column density along the line of sight, 
ctt is the Thomson cross-section, Lii nc is the line luminosity, and 
R is the typical size. (We defined Q = 10 X Q X in cgs units.) The 
energy-dependent opacity r 77 (E, Eu ne ) is obtained by multiplying 
tt by the angle-averaged 77 cross-section in units of err, which 
has a maximum of ~0.2 at three times the threshold energy (see 
I Aharonianl 120041 : PS 10). The opacities produced by two "lines", 
thc and th, were free parameters. We see that with the power-law 
null-hypothesis all spectra have very significant breaks and require 
non-zero absorption opacities (column 4 in Table [T). However, the 
residual value of \ 2 is still high for many spectra, especially for the 
average spectrum corresponding to the entire observational period, 
mostly because of the spectral curvature below the break. 

The simplest way to account for this curvature is to replace 
the power-law null-hypothesis by the lognormal distribution (or 
"log-parabola", whi ch is essentially the s a me function in a differen t 
representation, see Massaro et al. 2004|> lAckermann et al.l 1201 Oh . 
which we use in the decimal form 



EF E oc 10 



■log^iS/Bp^)/^ 



(3) 



Gamma-ray emission of 3C 454.3 3 



Table 1. Spectral properties of 3C 454.3. 



Interval 


Dates™ 
MJD 


Power-law 


x 2 

PL + DA 6 


/dot' 

Lognormal 


Logn + DA C 






MeV 






X 2 /dof^ 


A 
/A 


J^tOoH — J'Hr 1 Dy 




7 7/1 n 


37/1 1 


7 7/Q 




_.o 




7 + 1.0 
J - z -0.9 


7 ,+29 
zl -14 


10/ 1U 


B 


54759-55034 


34/12 


8.4/10 


16/11 


7.3/9 


■X 1+0.9 


0.68 


11±3 


4 5+ 13 


4+7.0 


7.5/10 


C 


55034-55159 


65/12 


6.8/10 


16/11 


3.7/9 


3.2 ±0.6 


2.4 


15±3 


4.0±0.9 


o ,+8.3 


3.8/10 


D 


55159-55199 


81/12 


6.7/10 


14/11 


2.9/9 


3 1 +1 '° 
J - 1 -0.4 


9.1 


49±6 


2.8±0.7 


4 6+30 


3.2/10 


E 


55199-55284 


63/12 


40/10 


4.4/1 1 


4.4/9 


1.6+.0.1 


3.0 


23±4 


4.0±1.0 


T. q+5.5 


19/10 


F 


55284-55344 


112/12 


23/10 


27/11 


17/9 


7 q+0.5 


7.6 


33±5 


2.6±0.6 




17/10 


G 


55344-55517 


85/12 


13/10 


18/11 


7.9/9 


i 7+0.7 


3.4 


19±3 


2.1±0.6 


3.it 2 2 : 


8.8/10 


H 


55517-55522 


79/12 


18/10 


11/11 


5.9/9 


7 fi +0.2 
zo -0.3 


55 


126±17 


1.6+.0.6 


7 ,+2.8 


6.8/10 


I 


55522-55587 


242/12 


33/10 


28/11 


10/9 


2.5+.0.2 


15 


60±5 


2.7±0.4 


c 7+2.6 
1.5 


11/10 


J 


55587-55652 


43/12 


11/10 


13/11 


7.2/9 


7 4+0.6 
z -^-0.4 


2.6 


26±6 


3 1+ 1 ' 3 
3 - 1 -l.l 


>8» 


7.3/10 


A+C+E+G 




315/12 


33/10 


42/11 


9.0/9 


2.7±0.2 


2.9 


20±2 


3.2±0.4 


fi2 +2. 2 
az -1.8 


9.0/10 


D+F+I 




565/12 


56/10 


71/11 


23/9 


7 O+0.3 
Z - S -0.2 


11 


48±3 


2.7±0.3 


7.2±1.5 


23/10 


A-J 


54684-55652 


1099/12 


85/10 


129/1 1 


23/9 


2.6±0.2 


3.9 


35±2 


2.7±0.2 


8.lii;§ 


23/10 



a The start and the end of the interval. "Power-law plus double absorber (DA) model. c Lognormal plus DA model. d Best-fitting values of <ri n for the 
lognormal + DA model. e Flux EFe at 1 GeV in units of 10~ 10 erg cm -2 s _1 . Best-fitting values of i? pca k, THe and tjj and the corresponding x 2 /dof for 
the lognormal (with cr\ n = 2.7) + DA model. 3 1<t lower limit. 
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Figure 2. (a) Spectral energy distribution (SED) of 3C 454.3 for ten intervals A-J. (b) SED of 3C 454.3 for intervals B and H of the lowest and highest fluxes 
as well as averaged over the whole observation period of 2.5 years (sum of intervals A-J). The solid curves show the best-fitting lognormal (<ri n = 2.7) + 
DA model. The dashed curve represents the fit with the lognormal distribution absorbed by BLR at a fixed ionization parameter log § = 2.5 (see PS 10). The 
dotted curves show the unabsorbed lognormal distributions. 



The resulting fits are shown in columns 5 (without absorption) and 
6 (with DA) of Table Q] We see a significant decrease in \ 2 com- 
pared to the power-law-based models. The DA model gives a sub- 
stantial reduction in \ 2 in all cases (except for interval E). These 
fits are superior compared to any simple phenomenological model 
considered earlier (e.g. lognor mal, broken powe r-law, or an expo- 
nentially cutoff power-law, see lAbdo et alj|201lh . 



If the spectral break is indeed caused by the 77-absorption on 
the He II LyC, its position should be stable. In that case we can 
expect that the break should become more significant in the aver- 
age spectrum corresponding to the whole observation period. In- 
deed, this is the case, the lognormal distribution without absorption 
fits the data (see Table [TJ at x 2 = 129 for 11 degrees of free- 
dom (dof). The DA improves the fit (see Fig.|2jb)), but the residual 
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Figure 3. Dependence of the width of the best-fitting lognormal + DA 
model on flux at 1 GeV. The parameters for the intervals A+C+E+G, A— 
J and D+F+I are shown by open circles. 



value x 2 /dof = 22.8/9 is still high. The absorption model for 
the BLR at a fixed ionization parameter log£ = 2.5 (see PS 10) 
gives xVdof = 27/10, with the quality of the fit being better 
(X 2 /dof = 7/7) in the range <15 GeV. The main contribution 
to the residual \ 2 comes actually from the high-energy tail of the 
spectrum and might have a simple interpretation: the absorber op- 
tical depth varies during the observation and the sum of the spectra 
with low and high opacities differs from the spectrum at a fixed 
intermediate opacity. 

The best-fitting value of a\ n is nearly constant ~ 2.5-3 (see 
column 7 of Table[T]and Fig. [3j, with the interval E being again an 
outlier. In order to avoid extra statistical noise, we fix eri n = 2.7 
and then probe the correlation of other parameters with the source 
brightness (see Table[T). We find a highly significant correlation of 
the peak energy Speak with flux, which is well represented by a 
power-law (Fig. |4j: 



Speak = (46 ± 2) F 



fl.60i0.04 



MeV, 



(4) 



where F = EFe at 1 G eV. This confirms a pre viously known 
hardness-flux correlation ( Acke rmann et al.l l2010b - We also ob- 
serve a marginally significant (2.5a) anti-correlation between the 
77-opacity by the He II LyC and the flux (see Fig. [5}: 

tho = (2.61 ± 0.22) - (1.15 ± 0.47) log F_ 9 . (5) 

The opacity th in the H I LyC does not vary significantly. 



4 DISCUSSION AND SUMMARY 

The BLRs around quasars emit a number of lines associated with 
very different ionization stages, with th e high-ionization lines be- 
ing broader dPeterson & Wander] Il999h . Reverberation mapping 
demonstrated a strong anti-correlation between the line width and 
the time delays of the lin e response to the continuum variations 
dPeters on & Wandel 2000). All this implies strong stratification of 
the BLR tha t extends over two orders o f magnitude in distance 
dKroiikll999llOsterbrock & Ferlandl2006l) . The BLR size, as mea- 
sured from the C IV A1549 line delays, scales with the accretion 
luminosity as Ren ~ 0.2L\j 2 pc dKaspi et"al] 2007). However, 
the delays in the hig h-ionization Hell A 1640 line are 3-5 time s 
smaller in Seyferts dKorista et alii 1 9951 ; IPeterson & Wandelll 1999b . 
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Figure 4. Dependence of the peak (in EFe) of the best-fitting lognormal 
(with (Ti n = 2.7) + DA model on flux at 1 GeV. The dashed line shows the 
best-fitting relation for intervals A-J given by Equation (|4j. The parameters 
for the intervals A+C+E+G, A-J and D+F+I are shown by open circles. 
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Figure 5. Dependence of the Thomson optical depth in the He II LyC (de- 
fined by Equation 0) on flux at 1 GeV. The model and notations are the 
same as in Fig. [4] The dashed line shows the best-fitting relation for inter- 
vals A-J given by Equation {5). 



while the Balmer lines give sizes 2-3 times larger. As the opacity 
in the GeV range depends on the line compactness (L/R), mod- 
els assuming that all l ines are produced at the same distance from 
the central source (e.g. lLiu & Baj2006l;lReImerl2007l) . strongly un- 
derestimate the 77-opacity by the relatively weak, high-ionization 
lines produced in a small volume (if the 7-rays are produced closer 
to the centre) and overestimate the opacity for multi-GeV photons 
from the strong, low-ionization lines which are produced much fur- 
ther away than the fiducial distance obtained from C IV. 

A high average 7-ray luminosity of 3C 454.3 du ring the last 
two y ears with flares up to L 1 ~ 2 x 10 50 erg s" 1 dAbdo et aU 
l201ll) provided an order of magnitude incr ease in photon stat istics 
compared to the earlier data analysed by Ab do et al. 1 d2009b and 
PS 10. In became even more evident that the spectral breaks at ~3 
GeV are real. The spectra at all fluxes are well described by a log- 
normal distribution absorbed by the He II and H I LyC at 54.4 and 
13.6 eV. The fits with this model are superior compared to the pre- 
viously used simple phenomenological models, especially at high 
fluxes with high photon statistics. 
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The BLR size corresponding to the He II and H I LyC can be 
estimated from Equation ([2}: 

R Hc « 0.22 L Hc ,44(W3) pc, Rh ~ 3.7 L HA5 ( m /7) pc. (6) 

IWills et all dl995h report the luminosity of 3C 454.3 in Lya of 
~ 10 45 erg s" 1 and in Hell A1640 of 6 x 10 43 erg s~\ As- 
suming equal luminosities in LyC and Lya, the resulting Rh is 
much larger than the usually quoted sub-pc size of BLR, but is 
consistent w ith the absence of Lya variability in powerful quasars 
dKaspi etalj|2007h . The Hell LyC luminosity is probably larger 
than 10 44 erg s _1 , but is model dependent and therefore the esti- 
mate of Rhc is less certain. A rather large ratio thc/th (between 
about 1/4 and 1) and the observed thc-Aux anti-correlation indi- 
cate that the 7-ray emission region lies close to the boundary of the 
He III zone at R < Rn B and moves out at higher luminosity. Ab- 
sorption troughs at 0.3-0.7 GeV due to the high-ionization C and 
O lines are potentially detectable, but require high photon statistics 
and better calibration of LAT. 

We note that the typical values of th c ~ 2-4 imply that the 
maximum optical depth due to the He II LyC photons (at about 7 
GeV - three times the threshold value of 4.8/(1 + z) GeV) varies 
around 0.4-0.8. The absorption by He II LyC is very significant and 
accurately describes the data in the 3-15 GeV range. At high en- 
ergies, the opacity decreases roughly inversely with the photon en- 
ergy and becomes negligible above ~300 GeV. The opacity there 
is dominated by the optical/IR lines which are mostly produced in 
the low-ionization zone of BLR at large distances. This strongly re- 
duces the opacity and can make the BLR transparent for the multi- 
GeV to TeV photons. 

A probable excess of photons above 15 GeV (especially for 
the average spectrum) over the prediction of the model with the 
best-fitting absorption indicates that the spectrum is a superpo- 
sition of emission states with different opacities due to the ex- 
tended or moving emission zone. This proposal is supported by the 
spectral variability seen during the latest flare in November 2010 
(Abdo et al. 2011), where photons above 10 GeV mostly arrive at 
the end of the flare, clearly indicating the position of the 7-ray emit- 
ting region further away from the black hole. 

Our findings that the spectral peak of the underlying contin- 
uum strongly correlates with the flux, with the width of the lognor- 
mal distribution being roughly constant, indicate a more efficient 
acceleration at high fluxes. These variations might not be consis- 
tent with the interpretation of the spectral curvature as the Klein- 
Nishina reduction of the Compto n scattering efficiency of the Lya 
photons dAckermann et al.l2010n . This in turn would imply that the 
seed photon energy density is not dominated by the hydrogen Ly- 
man continuum and line photons (in the high-ionization zone of the 
BLR the metal lines and the He II Lyman photons are actually more 
energetically important). Alternatively, the steady electron distribu- 
tion is modified by inefficient cooling in the Kle in-Nishina regime 
dZdziarski & Krolikll993l ; lModerski et alj2005h . 

The GeV breaks observed in bright blazars still do not answer 
the question what is the 7-ray production mechanism. It can be 
Comp ton scattering of the accretion disc dDermer & SchlickeiseJ 
Il993h or BLR photons dSikora et all 19941) . or even the synchrotron 
radiation of relativistic pairs. The efficient means for particle ac- 
celeration can be provided by the photon breeding mechanism 
dStern & Poutanenll2006ll2008() . 

A small distance to the beginning of the 7-ray emitting re- 
gion is supported by the fast variability obs erved from 3C 454.3 
dTavecchio et alj 2010; Bo nnoli et alj 1 20 lib . Association of the 
beamed 7-ray emission with a relativistic jet implies that the jet 



is already accelerated at a sub-parsec distance from the black hole 
(i.e. at thousand Sch warzschild radii from the ~ 1O 9 M0 black hole, 
iBonnoli et al.ll20lH) . This, in turn, argues against the jet laun ching 
mechanism by the accretion disc dBlandford & P ayne 1982), with 
the only known alternative being the extracti on of the black hole 
spin energy by the Blandford- Znajek process (Blandfo rd & Znaiekl 
1977; Ko missarov etafll2007l) . 
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